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This is an analytic study of the two-dimensional XY spin glass with ± J disorder. The Hamiltonian 
has a continuous spin symmetry and a discrete chiral symmetry, and therefore possesses, potentially, 
two different order parameters and correlation lengths. The cost of breaking the symmetries is 
probed by comparing the ground state energy under periodic (P) boundary conditions with the 
one under antiperiodic (AP) and under reflecting (R) boundary conditions. Two energy differences 
("domain wall energies") appear, AE AP and AE R , whose scaling behavior with system size is 
nontrivially related to the correlation length exponents. 

For a specific distribution of the ± J disorder we show that the chiral and spin correlation lengths 
diverge with the same exponent as T J. 0. The common exponent has a common cause, viz. the low 
reversal energy of domains of chiral variables. For general disorder we give a heuristic argumentation 
in terms of droplet excitations that leads to spin ordering on a longer, or equal, scale than chiral 
ordering. These results are in contrast with interpretations of Monte Carlo simulations. 

I. INTRODUCTION 

Due to the rotational symmetry of the XY model, any of its ground states is necessarily part of a continuum of 
ground states related by global spin rotations. It was first pointed out by Villain [1] that the ground state of an XY 
spin glass with random ±J interactions has, in addition, a twofold degeneracy. The two continua can be deduced 
from one another by a global spin reflection with respect to an arbitrary axis. They are characterized by opposite 
"chiralities" , that is, by an opposite sense of rotation of the spins as one moves around a plaquette of the lattice. 
At finite temperature one may have domains belonging to different ground states ( "chiral excitations" ) , so that each 
plaquette has to be characterized by its own chiral variable. 

A number of recent theoretical papers have focused on the role of these chiralities [2-6] . The motivation for this 
interest is the question what the lower critical dimension dg is of the XY spin glass, and by which mechanism this 
system orders just above dg. Ozcki and Nishimori [7-9] have developed arguments for the bound dg > 4 (see also 
Schwartz and Young [10]), which is supported by substantial numerical evidence (see references in [9]). However, all 
these authors consider the conventional Edwards-Anderson order parameter associated with the XY spins. Villain's 
[1] discovery of a discrete symmetry therefore naturally led to the idea, due to Kawamura and Tanemura [2, 3], that 
in fact dg might be less than four. We recall briefly how this idea has been investigated in recent years. 

One can define two correlation lengths £ c and £ s , associated with the chiral variables qr (where r denotes the center 
of a plaquette) and the spin variables (where i denotes a lattice site), respectively : 

< Qrq r+R > 2 ~ e- R/ t° and < S t ■ S i+R > 2 ~ e~ R/i ' {R -► oo) . (1.1) 

Here < • • • > denotes the thermal and 7TT the disorder average. For the low temperature behavior of these correlation 
lengths one expects for d < dg : 

MT)~T-"< and &(T) ~ T~ v ' (T — > 0) . (1.2) 

One method to investigate the relation between £ c and £ s is to determine the correlation length exponents v c and 
v s via finite-size scaling of "domain wall" energies. In simpler cases where there is a single symmetry and a single 
correlation length with exponent v, this is done by finding the ground state energies for two appropriately chosen 
distin ct boundary conditions. The energy difference, AE, will then generally scale with the linear system size N as 
\AE\ ~ N- 1 '". 

In the present case, with both a continuous and a discrete symmetry, the periodic (P) "reference" boundary 
conditions can be changed in two ways : to antiperiodic (AP) and to reflecting (R) ones, this latter possibility having 
been first suggested by Kawamura and Tanemura [3] . Hence there are two energy differences to be considered, AE AP 
and AE R . Across an antiferromagnetic seam two spins interact with each other's image under reflection in the origin 
(of spin space). Hence AP boundary conditions probe the spatial rigidity of spin order, and the spin correlation length 
exponent v s is given by the usual relation 

\AE AP \ ~ N- 1/v ' (1.3) 

Across a reflecting seam two spins interact with each other's image under reflection about a fixed axis (in spin space). 
Since the chiral variables appear only when the "spin waves" are integrated out - and notwithstanding the obvious link 
between R boundary conditions and chirality - it is not a priori clear exactly how these boundary conditions probe 



1 



the chiral order. In fact the relation |A£^| ~ N^ 1 / 1 ' that one might naively expect is certainly not correct. The 
reason is that R boundary conditions, upon closer inspection, appear to affect also the spin wave degrees of freedom. 
This was first realized with remarkable intuition by Kawamura and Tanemura [3] , who then proceeded and extracted 
the chiral correlation length exponent v c from a heuristic expression involving both AE AP and AE R . 

On the basis of Monte Carlo determinations of the ground states of two- and three-dimensional systems these 
authors conclude that in two dimensions v s = 1.2 ± 0.15 and v c = 2.6 ± 0.3. This implies that the chiral variables 
have their own length scale and order in larger domains than the spin variables. Similar conclusions, namely v s m 1 
and v c ss 2, were reached by Ray and Moore [5], also on the basis of Monte Carlo simulations and a finite size scaling 
analysis. These studies suggest that as the dimension d is increased, v~ may vanish before v~ does, so that just 
above dt there would be a phase with long range chiral order but exponentially decaying spin-spin correlations. Indeed 
this is precisely the scenario that Kawamura and Tanemura [2,3] find in Monte Carlo simulations of the two- and 
three-dimensional XY spin glass. In recent work on the two- and three-dimensional Heiscnbcrg spin glass Kawamura 
[11] arrives at fully analogous conclusions. If these are accepted, the XY and Hcisenberg vector spin glasses have their 
lower critical dimensionality between d = 2 and d = 3. 

In order to achieve a better understanding of the interplay between chiral and spin variables, Ney-Nifle, Hilhorst, 
and Moore [6] recently performed an analytic study of the random ± J XY spin glass (in its Villain [12] formulation) 
on the one-dimensional ladder lattice (see also Morris et al. [13]). Their conclusion is that the chiral correlation length 
£ c and spin correlation length £ s diverge for T j with the same exponent v c = v s (of which they determine the value 
1.8999 ...). The purpose of this work is to extend the methods of Rcf. [6] to d = 2 and to confront our results with 
those of Rcfs. [3] and [5]. 

Our starting point in Section 2 is the Villain [12] formulation of the XY model with random ± J interactions. This 
formulation has the advantage that the chiral variables can easily be defined and can be decoupled from the spin 
waves. We transform the XY Hamiltonian to a Coulomb gas Hamiltonian with charges qr that play the role of the 
chiral variables. The transformation has been known [14-16] since long, including the fact [1, 17] that the charges take 
half-integer values on the frustrated plaquettes and integer values on the others. Here we take exactly into account 
all finite size effects, essential for the finite size scaling analysis that follows. Our result for the Coulomb gas partition 
function in the case of periodic (P) boundary conditions, and for an arbitrary realization of the disorder, is given by 
Eqs. (2.38). The effective Hamiltonian contains, in addition to the Coulomb interaction, a coupling term between the 
total electric dipolc moment and the boundary conditions. This additional term was obtained via a different approach 
in a very recent paper by Vallat and Beck [18], and its one-dimensional equivalent had appeared in [6]. In Sections 
3 and 4 we go further and discuss how these equations are modified in the case of antiperiodic (AP) and reflecting 
(R) boundary conditions, respectively. In the limit T — > the equations reduce to expressions for the ground state 
energy, including all its finite size corrections, provided the ground state itself is known or can be plausibly guessed. 

Since the ground state problem for an arbitrary realization of the disorder cannot be solved, we treat in Section 
5 a more restricted two-parameter subset of realizations, in which the frustrated plaquettes are placed in a random 
rectangular array with infinite range correlation along the y direction. (Disorder of the unidircctionally infinite-ranged 
type was also the first one considered for the random bond Ising model [19].) We find a regime of main interest in 
parameter space in which the XY spin glass has a zero-temperature transition with 

v s = v c (1.4) 

There is a second regime, with a strongly anisotropic spatial distribution of the frustrated plaquettes, in which the 
model has a low temperature phase with long range chiral correlation and, very plausibly, power law decay of spin 
correlations. This is behavior analogous to that of fully frustrated two-dimensional XY models [20-22], on which there 
exists a large literature. 

The result (1.4) in the first regime is different from the scenario proposed by Kawamura and Tanemura [2, 3] 
and by Ray and Moore [5]. The basic mechanism responsible for (1.4) is that the ground state can accommodate 
to a changeover from P to AP boundary conditions by means of the formation of a chiral domain wall, which is 
energetically lower lying than the continuous spin wave deformation that naturally comes to mind. On the basis of 
this example alone we cannot rule out the possibility that our conclusion, Eq. (1.4), is valid only within the restricted 
class of disorder realizations. 

In order to treat the case of general disorder, we construct in Section 6 a heuristic theory based on the same 
mechanism. It involves the lowest lying excitations of the standard Coulomb Hamiltonian which are assumed to be 
collective charges reversals. We find that 

v s > v c . (1.5) 

This is compatible with (1.4) but contradicts the results of Refs. [2,3,5] where the opposite inequality holds. Further- 
more, we show that if Eq. (1.5) would holds as a strict inequality, one cannot extract the exponent v c from AE R and 
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AE AP . To find v c , an appropriate quantity would then be the energy difference between the two ground states of the 
Coulomb Hamiltonian with periodic and reflecting boundary conditions, without any additional terms. 

Our conclusion is that there is no evidence that chiral order extends on a longer length scale than spin order. 



II. PARTITION FUNCTION WITH PERIODIC BOUNDARY CONDITIONS 

A. The Villain XY model with ± J interactions 

We consider an XY model on a finite square lattice, periodic in both directions, with sites i = (m, n), where 
m = 1, ■ • • , M and n = 1, • • • , N. Each site i is occupied by a two-component unit vector or "spin" whose angle 
(f)^ with a reference axis takes values in (— n, n]. In the ferromagnetic Villain model [12] two spins (f>^ and <f>j linked 
by a nearest-neighbor bond <i,j> have the Boltzmann weight 

e J = > e 1 J . (2.1) 

n— — oo 

Here J > sets the energy scale. The relation (2.1) implies that the Villain interaction V{<p) depends on (3J ; when 
(3 becomes large, only one term on the RHS of (2.1) will dominate, V(</>) will tend to 

V((j>) = 4? for |0| < 7T, (3 = oo (2.2) 

and we may interpret (3 as the inverse temperature l/k B T. For all (3, the sum on n in (2.1) guarantees the periodicity 
property V(ip) = V((p + 2ir) and the symmetry property V(ip) = V(— ip). The antiferromagnetic Villain model has 
V(<fii — 4>j — 7r) instead of V(<p^ — <j>j). 

Here we wish to consider an XY spin glass with randomly fcrro- or antiferromagnetic Villain interactions, that is, 
with the Hamiltonian 

H = J J2 V^i-^j-Hj) ( 2 - 3 ) 
<i,j> 

where the sum runs over all nearest-neighbor bonds of the periodic lattice, and the TT^j are quenched random variables 
such that 
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' with probability ^ (2.4) 



% 3 Kf ^—J 2 

The expression for the canonical partition function of this XY spin glass then is 

z m .n= / Q<% e <l,3> ( 2 - 5 ) 

in which the n^j may be seen as additional dynamical variables, and the argument of the exponential as a new effective 
Hamiltonian. We shall henceforth write 

u ij =n ij /(27r) , (2.6) 
which is integer (half-integer) when the interaction between ip^ and ipj is ferromagnetic (antiferromagnetic). 

B. Transformation to a Coulomb gas 

1. From XY spin glass to a SOS model 

The transformation from an XY model to an SOS model is well-known [14,15] and has been extended to various 
types of random XY models. Here we carefully study the finite size effects, that determine the ground state energy 
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differences under different boundary conditions. In (2.5) we wish to carry out the integrations on the MN variables 
</9j. To that end we arbitrarily select a lattice site io and transform to the new variables of integration 



(2.7a) 



(2.7b) 



These are 2MN + 1 in number, and there exist MN + 1 relations between them. In order to formulate these, let r, s, 
• • • be the position vectors of the plaquette centers. We adopt the convention that in <i,j > the site j is to the right 
of i (for a horizontal bond) or above i (for a vertical bond), which naturally extends across the periodic boundaries. 
We define furthermore a sum on the bonds < > surrounding a plaquette r by the diagram of Fig. 1 



V43 



¥>14 




<P23 



1 ¥12 2 

Figure 1: Diagram indicating the sign convention in the plaquette sum of Eq. (2.8). The indices 1, ■ • • ,4 are shorthand for ii, • • • , 14. 

together with the formula 



^ <Pa' = ¥>12 + ^23 - </?43 - ifilA 



The hat on the summation sign is a reminder of the sign convention in the RHS of (2.£ 
then satisfy the relations 



(2.8) 



). The variables ip^j of (2.7b) 



tp^t = mod 2tt for all r 



(2.9a) 



These are MN relations of which only MN — 1 are independent. The two remaining relations correspond to loops 
around the torus and are 



M 



X! Vij = X! P(m,m),(m+l,ni) = mod 2lT 



(2.9b) 



m— 1 
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N 



^2 V(mi,n),(mi,n+1) = mod 2?r 



(2.9c) 



n=l 



where ii = (m\,n\) is another arbitrarily selected lattice site. Upon introducing the variables of integration (2.7) in 
Eq. (2.5) we must represent the conditions (2.9) by delta functions and find 
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Z M M = 2* JZoU<i,j> dl Pij Ur^ro 6 (£ <Pij mod 2n 

xS (j2 X( Pij mod 27r ) s (j>2 y mod 27r 
-PJ £ ( 



(2. 10) 



x cxp 



<ij> 



in which To is an arbitrarily selected plaquette, and where a trivial factor 2ir comes from the integration on ipo. The 
delta functions may be represented by a Fourier sum with the aid of 



S (x mod 2tt) = ^ 6{x - 2tt I) = (27T)" 1 ^ 



(2.11) 



which requires the introduction of a set {rir} of summation variables for the plaquettes, and of two summation 
variables n x and n y for the loops around the torus. One finds 



/oo 
II d nj exp 
-oo . . 



{ny} n x n y 



<ij> 
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x exp \inJY^ ¥>ij + in v Y^ \%j 



x exp 



<*.j> 



(2.12) 



where the prime on the summation sign indicates the restriction to n Ta = 0. The integrations on the ip^j can now be 
performed. 

Let the geometric relation between the pair of lattice sites < i,j > and the pair of plaquette centers < r, s > be 
as in Fig. 2. 



i •- 



-• 3 



r 9- 



-• s 



Figure 2: Geometric relation between the pair of lattice sites < i,j > and the pair of plaquette centers < r, s >. 

It is furthermore useful to define 



x(y) 

T rs 



| 1 if <i,j> is part of the loop (2.9b) (the loop (2.9c)) 
1 1 



otherwise 



(2.13) 



Henceforth we shall generally write n r s instead of n^j when < i, j > and < r,s > are related as in Fig. 2. With this 
notation the result of the ip^j integrations in (2.12) is 
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z M , N = (2pj)- MN j2T,H ex p 



n y {n r } 



* E 



> ttts I n r - n s + 

<r,s> \ 



7 y T rs n a 



x exp 



<r,s> 



rir — ns 



E 



(2.14) 



Eq. (2.14) represents the partition function of a solid-on-solid model (or "column model" ) with Gaussian interaction, in 
which, due to the randomness in the original XY model, some terms in occur with negative sign. In the nonrandom 
case with all ttvs — 0, we recover the partition function of the well-known discrete Gaussian model, summed on 
different step boundary conditions represented by the variables n x and n y . 

2. From SOS model to Coulomb gas 

The transformation from an SOS model to a Coulomb gas is also well-known [16]. We have to apply it here to the 
random SOS model of Eq. (2.14), taking properly into account again all finite size effects. For a function f(n) of an 
integer variable n one has 



E /(»)= E / dve^f{v 



(2.15) 



Applying this identity to the sum on the {n r } in (2.14) we get 



(2(3 J) 



-MN 



/Seee fu 



dv r x 



n x n y { qj ,} 



x exp 



2ni ^2 1r v r + 2m\v + i £ n r s I vr - v s + £ t£ s n a ] 
r <r,s> \ a ) 



<r,s> 



- ^ + E r ^ n °^j 



(2.16) 

in which = vr a and the integral on A takes care of the condition rir = in (2.14). The integrations on the 
are Gaussian. The most convenient way to carry them out requires some preliminaries. We introduce the coordinate 
representation r = (x, y) with x = 1, • • • , M and y = 1, • • • , N, and the periodicity condition that (x + M, y) and 
(x, y + N) also denote the plaquette center r. For horizontal and vertical pairs < r, s > we shall have, by convention, 
s = (x + 1, y) and s = (x, y + 1), respectively. For each plaquette center r we define a frustration variable p r by 



2tT p r = £ 7T = 7T 01 - 7T 2 - 7T 30 + 7T 4 

where the notation is as in Eq. (2.8), together with Fig. 3. 



(2.17) 
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Figure 3: Diagram indicating the sign convention in the definition (2.17) of the frustration variable pr- The indices 0, 1, 
for the plaquette centers r, ri, ■ ■ ■ , 1-4. 



, 4 are shorthand 



Let furthermore 



(2.18) 



<r,s> 



where the subscript a — x (or a = y) denotes restriction to vertical (or horizontal) bonds (recall that a horizontal 
pair < i,j > corresponds to a vertical pair < r, s > and vice versa). At this point it is useful to fix the positions of 
the two loops around the torus by the specific choice 



T rs = S v ,n for <r,s >=< (x, y), (x, y + 1) > 



T rs = S x,m for < r,s >=< (x,y), (x + l,y) > 



We pass in Eq. (2.16) from the v r to new variables of integration v' r defined by 



v' r = vr + M 1 xn y + N l yn x 



(2.19) 



(2.20) 



where as before r runs through (x, y) with x = 1, • • • , M and y = 1, • • • , N. After rearranging terms in the exponential 
one finds for the partition function (2.16) 



Z MtN = {2f3J)- MN / dX exp [- 27T i\(n x + n y ) + i (M^UyHy + A^ _1 II x n x ) - (4/3 J)" 1 (N^Mnl + M~ 1 Niiy)] 



2 n x n 



x Y ex P 



yn x 



x Z' 



M,N 



with 



/oo 
\\dv' r exp 
-OO T> 



<r,s> 



We shall abbreviate 



(2.21a) 



(2.21b) 



Q r = q r +p r + \S rro . 



(2.22) 
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The integrals on the v' r in Eq. (2.21b) are now easily carried out with the aid of the Fourier variables 



(2.23) 



and analogously defined Q^, where 



k = (k x ,ky) = 2n(M- 1 K X7 N^ 1 K y ) , 



(2.24a) 



k x = 0,1,---,M-1 , k v = 0, N- 1 . 



(2.24b) 



It is useful to make the specific choice ro = (M, N). Taking properly care of the integration on Pq, which is exceptional 
and gives a factor 5(Qq), one finds from (2.21b) 



with 



z 'm,n = {MN)*5 ir + A J C m ,n exp 



fc#0 



A fc = 4 | sin 2 ^ + sm 



2 n-j/ 



(2.25) 



(2.26) 



c m ,n = n (*x0 j /*kY 



(2.27) 



Upon inserting (2.25) in (2.21a) one may carry out the A integration with the result 

Z m ,n = {2f3J)- MN {MN)ic M ,Nj2J2 e *V [* ( M ~ lll y n y + W^n*) - (4/3J)" 1 (N^Mnl + M" 1 ^ 2 )] 



x ' ex P - 2lTi E ( M ~ lxn v + ^~ V*) - (3H C ({qr + Pr}) 
{ qr } L r 

Here the prime restricts the summation to neutral "charge" configurations {qr}, i-e. satisfying 

r 

and He is the Coulomb Hamiltonian 

H c ({qr + Pr}) = 8ir 2 J ^ E U m ,n(t - r ')(q r + p r )(qr> + Pr') 

r r ' 

where Um n is the Coulomb potential on a periodic lattice, 



(2.28) 



e - t k.r _ j 



fc#0 



(2.29) 



(2.30) 



(2.31) 
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This function has the periodicity properties 



Um,n(t) = U m ,n (r + Mei) = U m ,n (r + Ne 2 ) (2.32) 

where e\ = (1,0) and e 2 = (0, 1). 

The last step needed is to transform the sums on n x and n y in (2.28) according to Eq. (2.15), by introducing 
continuous variables v x and u y , and new summation variables q x and q y . The integrations on v x and v y are again 
Gaussian and easily carried out. The result can be slightly rewritten with the aid of the relations 

M~ x (^l y - 27r^x gr ^ = -2-kM~ x P x + 7r y (2.33a) 



-27r^j/g r ^ =-2^-^+^ 



(2.33b) 



where ir y and 7r x are sums along loops around the torus, 



N M 

n V = ^2 n (M, y ),(l,y) , 7I"x = ^7T( XiJV ),(x,l) (2-34) 

y—1 x—1 

and P is the electric dipole moment 

P = (P x ,P y ) = J2 r (<lr+Pr) ■ (2.35) 
r 

The final result for Zm,n then becomes 

Z m ,n = (2(3J)- mn (MN)^C m ,n ]T 'Qn/m (N^Py - n x /(2ir)) Q M/N {M~ 1 P X - ir v /(2n)) e~ mc (2.36) 

where 



e a (u) = (47r/3Ja)3 ^ exp [-4ir 2 f3Ja(q - u) 2 ] (2.37) 

g= — oo 

has the periodicity property Q a (u) = o (u + 1), and the prime refers to the charge neutrality condition (2.29). 

C. Summary 

In view of the relatively technical character of the above transformation, and for easy later reference, we summarize 
here the result. The partition function Zm,n (see Eq. (2.5)) of the ±J Villain XY model (Eqs. (2.3) and (2.4)) can 
be expressed as the partition function of a system of Coulomb charges {qr} on the dual lattice, 

oo 

Z m ,n = 2k(2PJ)- mn+1 {MN)^C M m J2' e ~ 0n{9x ' 9y ' i9ri) ( 2 - 38a ) 

Qx,q v =— oo {q r } 

where 

M 

n{q x ,q y , {q r }) = 4tt 2 J7VM" 1 [q y + AT 1 ^ y(q r + Pr) + (27T)" 1 ^ 7r (XiJV ),( x ,i)] 2 

r x=i 
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N 

+4n 2 JN- 1 M [q x + M- 1 Y, x(q r + Pr) + i^Y 1 ^(M, y ),(i, y )] ' 

r y=l 

+87r 2 J^^f/ M> Jv(r - r')(q r +Pr)(q r > +P r >) (2.38b) 
r r > 

The Coulomb potential Um,n and the constant Cm,n are given by (2.31) and (2.27), respectively, together with 
(2.26) and (2.24). The qr run through all integer values subject to the charge neutrality condition (2.29). The sums 
on r and r' run through (x, y) with x = 1, . . . , M and y = 1, . . . , N. Each nearest neighbor bond (r, s) on the charge 
lattice is dual to a nearest neighbor bond of the original spin lattice, and the disorder variables Tr rs that appear 
in (2.38) are equal to the corresponding TTjj- The frustration variables pr are defined in terms of the ■nrs by the 
diagram of Fig. 3; they are half-integer for a frustrated plaquette, and integer for an unfrustrated plaquette. Even 
though the first two terms in (2.38b) seem to favor a definite coordinate representation of the lattice, one may verify 
with the aid of some algebra that Zm,n is invariant under translation of the origin, as of course it should be. 

The transformation from the ferromagnetic 2D XY model to the Coulomb gas Hamiltonian Tic, including the 
charge neutrality condition, has been known since long [15]. Similarly the replacement of the charges qr by qr + pr 
in the random case has been known since Villain [1] and was treated in a more general context in Ref. [17]. Eqs. 
(2.38) show how on a finite lattice the charge configurations receive a supplementary weight, represented by the first 
two terms in (2.38b), that depend on the electric dipole moment, in agreement with [18]. Whereas this weight plays 
no role for a bulk system in the thermodynamic limit, it is essential for finite size scaling considerations, and therefore 
(see Eqs. (1.2) and (1.3)) for the determination of the low temperature behavior of the correlation lengths. This fact 
will be exploited in the example of Section 5 and in the general heuristic theory of Section 6. 



III. ANTIPERIODIC BOUNDARY CONDITIONS 



In view of the remark preceding Eq. (2.3), it is very simple to change the periodic boundary conditions in the 
system discussed above into antiperiodic boundary conditions (we shall always take these along the seam joining the 
Mth and the 1st column). It amounts to changing n^j into ir^j + n on the antiferromagnetic seam. This just means 

drawing another member of the class of random systems under consideration. In particular, frustrated (unfrustrated) 
plaquettes remain frustrated (unfrustrated). The corresponding antiperiodic partition function Z^f N differs from Eq. 
(2.36) only in that tt x /(2it) is replaced by tt x /(2tt) + \. This difference is at the origin of the finite size effect of 
interest. 

The following example, although trivial, shows how ground state energy differences are extracted from the final 
equations of Section 2. We compare the ground states of an XY ferromagnet with periodic and antiperiodic boundary 
conditions. The energy difference AE AP is that of a spin wave of wavelength 2M, which can be written down 
immediately : 

AE AP = NMj(^-] = tt 2 ^J . (3.1) 
\2M J M K ' 

To see how the same result can be obtained from Eq. (2.36) we use that for both types of boundary conditions all pr 
vanish, and in the ground state all qr as well ; furthermore that the periodic boundary conditions have ir x = and 
the antiperiodic ones tx x ~\. This gives 

AE AP = - lim hogQ N /M(b (3-2) 

which in view of (2.37) exactly coincides with (3.1). 

The functions in Eq.(2.36), or equivalently, the first two terms in Eq.(2.38b), can be traced back mathematically 
to the global constraint on the spin variable differences fij- This example shows that they represent the energy 
of a continuous spin wave deformation forced into the system by the boundary conditons. We shall extend this 
interpretation to the case of general disorder, and speak of a "global spin wave" . 

Ground state energy differences between P and AP boundary conditions for more complicated situations can also 
be obtained from (2.36), or (2.38), at least in those cases where the ground state of the charge system is known or 
can be plausibly guessed. A nontrivial example is discussed in Section 5. 
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IV. REFLECTING BOUNDARY CONDITIONS 



Reflecting boundary conditions are introduced into the XY Hamiltonian by letting the spins <S(i )n ) in the first 
lattice column interact with the images of the S(M,n) in the last column under reflection about an axis in spin space. 
The counterpart of the periodic Hamiltonian (2.3) is 

W K ^E re9 v Vj Hj) +J Y. excy + <Pj Hj) ■ (4-1) 
<i.j> <i,j> 

Its partition function will be denoted by Zfy N , the upper index R indicating, here and henceforth, reflecting boundary 
conditions. The subscript "reg" in (4.1) refers to the regular terms and the subscript "exc" to the exceptional ones, 
modified by the reflecting boundary conditions. 



Transformation to a Coulomb gas 

The conversion of the Hamiltonian (4.1) into a Coulomb gas Hamiltonian proceeds via the same succession of 
transformations as in the case of periodic boundary conditions. However, several differences occur, and we shall 
indicate the main modifications below. 

(i) When transforming to the variables of Eq. (2.7), it is convenient to choose Zq = (1, 1). The expressions 
for tpj + ipj in the exceptional terms in (4.1) then become 

n 

<P(M,n) + <P(l,n) = <P(M,n),(l,n) ~ 2 1(1,1-1), (1,1) + 2 <P0 (4.2) 

t=2 



(n=l,2,...,A) . 

(ii) The integration on ipo can be carried out only after the introduction of the plaquette and loop variables 
{nr}, n x , andn y , and leads to the result 27r5„ xi o- The integrations on the <P(M,n),(i,n) an d <P(i,n-i),(i,n) are exceptional. 

(iii) The discrete Gaussian partition function, obtained as an intermediate result, now becomes 



7 R 



(2/3J)-^££'exp 

n y {n r } 



i ^ "rsinr - n a ) - (4/3J)" 1 ^ reg (n r - n s f 



L <r,s> 



<r,s> 



x cxp 



^ y exc 

<r,s> 



TTr,s(-n r -n s + n y ) - (4/3 J)" 



<r,s> 



'(nr - n s - n y y 



(4.3) 



in which, as before r = (x, y) with x — 1, • • • , M and y = 1, ■ • • , N, and the exceptional terms are the horizontal bonds 
linking the Mth to the first column ; we have chosen t^, s = 1 if < r, s > is exceptional (and hence t^, s = if < r, s > 
is regular) ; and the prime restricts the summation to configurations {nr} with n/M,N) — 0. The exceptional terms 
show that between the columns M and 1, in addition to the step boundary condition represented by the variable n y , 
a reflection is imposed with respect to the zero level of the column heights. 

(iv) Upon continuing the succession of transformations, one finds that the integral on v y gives 2S(2q y + 

+ A). After integration over A this becomes 2S(2q y — J^9r , 0) where S(a, b) = <5 a ^, and after summation on q y 
r r 
one gets 25{ S y^ j qr mod 2, 0), i.e. the sum of the charges qr should be even, 
r 

(v) The relevant wavevectors are now 



k = (k x , k y ) = 2n ^Al- 1 (^k x + , A^ 1 *^ 



(4.4) 
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with k x and n y as in (2.24b) ; this set of wavevectors will be referred to by the upper index R on summation and 
product signs. We define in particular 

CM,N = H R (^l3J/\ k ) i . (4.5) 
k 

(vi) The variable qr — S x .M^rs/^, which appears in the Coulomb Hamiltonian, may be renamed q r by a shift 
of the variables q(M,y)- 
The final result is 

Z^ N = 2(2(3 jy MN C^ N Par (4-6) 

{ir} 

where the subscript "par" refers to the parity condition 

(^2 1r + TTy/vj mod 2 = , (4.7) 

and where 

Wc (tir+Pr}) = 87T 2 jJ2J2 U MM r ~ r ') (<lr + Pr) {qr >+Pr>) (4-8) 

r r ' 

with r and r 1 in the range (x, y), x = 1, . . . , M, y = 1, . . . , N, and with 



2MN^ X k 



This function has the (anti-)periodicity properties 



U^ N {r) = U« iN (r + Ne 2 ) = -U^ N (r + Me,) . (4.10) 

Hence two charges at a fixed finite distance on opposite sides of the reflecting boundary interact with each other's 
charge conjugated image : reflecting boundary conditions for the XY spins lead to charge conjugating boundary 
conditions for the Coulomb charges. 

We now comment on this result and compare it to its counterpart for P and AP boundary conditions Eq. (2.36). 
The fact that in the partition function (4.6) no charge neutrality is imposed becomes understandable if one writes the 
interaction U R as 

Uf t , N {r) = U 2M , N {r) - U 2M A r + Me i) (4- 11 ) 

which is easily checked. It is as though we have a double system, of size 2M x N and with periodic boundary 
conditions, in which every charge at a site r of the original system is paired up with an image charge, equal but of 
opposite sign, at the corresponding site r + Me\. Using Eq. (4.11) and obvious symmetry properties one can make 
this explicit by writing the Hamiltonian (4.8) as 



n£({qp+Pp}) =4 7 r 2 J^^(7 2 M^(p-p')(9p+Pp)('7p'+Pp') , (4-12) 

P P ' 

where p and p ' run through the 2M x TV lattice and it is understood that corresponding sites carry opposite charges. 
The extra prefactor | that (4.12) has with respect to (4.8) indicates that the energy density has been smeared out 
over twice the volume. The self-interaction present in (4.8) has become the interaction energy between a charge and 
its image in (4.12). 
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V. AN EXAMPLE 



A. Specific expressions for the disorder. Ground state and ground state energy with P boundary conditions. 

As an example we consider an M x N lattice containing a rectangular array of M' x N' frustrated plaquettes, 
consisting of the points r a ,b = (x a , yb) with a = 1,2,..., M' and b = 1,2,..., N', where yb = bN/N' = b£ y with £ y 
an even integer, and the Xb are drawn randomly and independently with a density M'/M = l/£ x . The density of 
frustrated plaquettes is therefore l/£ x £ y . (Note that drawing the frustrated plaquettes randomly is not the same thing 
as drawing the negative bonds randomly [23, 24]; the latter procedure creates pairs of frustrated plaquettes, leading 
to a Coulomb gas with quenched dipolcs rather than quenched charges). We let Q% denote the ground state value of 
qr + Pr ■ The arguments of this section rest on the hypothesis that we can find the true ground state. For M' and N' 
even we may safely assume that the ground state has the chessboard array of charges 

«?...»> = 4 ( - i)o+6 (fu) 

and zero charges on all the other plaquettes. The ground state value of H c therefore is, from (5.1) and (2.30), 

H c ({Q° r }) = ^ 2 JJ2 E U M,N(ra,b - r a , M )(-l) a+a ' +b+b ' (5.2) 

a, b a' ,b' 

Upon using the explicit form of the lattice Coulomb potential, Eqs. (2.31) and (2.26), we can rewrite (5.2) as a sum 
of interaction energies between the charge carrying lattice columns, 

n c ({Q° r }) = ^J£- 2 NY,{-l) a+a ' E W*a - <) (5.3) 

a. a' v 

where v runs through the values ±5, ±|, . . . , ±{^£ y — 5), and in which, with c„ = cos2irv£ y 1 , 

= M S 4-2cos/c x -2c„ ^ 

* \ [(2 - c v f 1] - 1/2 {2 -c- [(2 - c v f 1] 1/2 }N 

the last step involving the limit M — > 00 at finite \x\. This expression takes the much simpler form 

V v (x) - (2Tr\u\)- 1 £ y exp(^2Tr\iyx\£- 1 ) (5.5) 

in the limit £ x ,£ y — > exo with £ x /£ y fixed and x of order In this limit, also, the sum on v in (5.3) can be evaluated 
exactly and one finds 

H c ({Q° r }) - -KJl- 1 NY,{-V a+a ' 'logtanh(7rK - ^K" 1 ) (5.6) 

a, a' 

This shows that the interaction energy is negative between two neighboring columns (a and a + 1), that it diverges 
logarithmically at small distances \x a — x a '\, as expected, and tends to zero as exp(— ir\x a — x a > l^ 1 ) for large \x a — x a >\, 
so that it is short-ranged. 

One remark is in place. The limit £ x , £ y — > 00 corresponds to a dilute array of frustrated plaquettes. These must be 
joined pairwise by long ladders of negative bonds - for which there is no physical reason. The analysis that follows 
will use Eq. (5.6); however qualitatively similar results can be obtained from Eqs. (5.3) and (5.4), where £ x and £ y are 
arbitrary, the key point being that the large distance behavior of the column-column interaction is still determined 
by the v — ±5 terms in Eq. (5.3). We shall come back to the case of general £ x and £ y in the discussion in Sec. 6. 

B. Comparison of P, R, and AP boundary conditions 

We now wish to consider the ground state energies and their differences for the three different boundary conditions, 
P, AP, and R. 

(i) Periodic boundary conditions 
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For the frustrated plaquettes located on the array defined in the beginning of this section, we can always choose 
the negative bonds such that n x = n y = 0. Since the ground state {Qr} has P x — P y = (see Eq. (2.35)), the first 
two terms in (2.38b) are minimized by q x = q y = and so do not contribute to the ground state energy, which is 
therefore entirely given by (5.6). 

One possible type of excitations from the ground state are those consisting of reversing all charges Q% in one column 
a. The excitation energy is of the order of 

AE(Z a ) ee 2^- 1 Anogtanh« Q ^ 1 ) (5.7) 

where £ a is the distance between x a and the closer one of ir a _i and x a+ \. Because of its exponential decay with 
\ x a — x a i\, the interaction between a and more distant columns does not change this estimate. 



(ii) Reflecting boundary conditions 

We first of all have to find the a priori unknown ground state. To this end we consider again the chessboard 
configuration of charges of Eq. (5.1) and calculate its energy H^({Qj.}) under reflecting boundary conditions, 
starting from Eqs. (4.8) and (4.9). The result is again Eq. (5.6), but with charge conjugating boundary conditions, 
the interaction terms between columns a and a' on different sides of the reflecting boundary have an extra minus sign. 
Since the number M' of columns was supposed even, this means that there is somewhere a mismatch consisting of two 
neighboring columns a and a + 1 that are identically instead of oppositely charged. The energy cost is ~ AE(£), with 
£ the distance between these columns, and can be minimized by shifting the mismatch to the largest interneighbor 
distance. The typical maximum distance that can be found is 

6nax-41og(^)+0(l) (M — > oo) (5.8) 

and therefore the ground state energy difference AE R between reflecting and periodic boundary conditions is of the 
order of Ai?(£ m ax), which leads to 

We can now link the finite size scaling of the ground state to the low temperature behavior of the correlation 
length as indicated in the introduction. Since the mechanism leading to Eq.(5.9) is the reversal of large domains of 
chiral variables, one deduces that in case nl x j£ y > 1, there is in the thermodynamic limit M ~ N — > oo no chiral 
order at any finite temperature T, and that for T [ the chiral correlation length £ c diverges with an exponent 
v c = 1/(1 — n£ x /£ y ). In case n£ x /£ v < 1 (which corresponds to relatively strong spatial anisotropy), one deduces 
similarly that the system has a low temperature phase with long range chiral order. We shall come back to both cases 
after discussing AP boundary conditions. 



(iii) Antiperiodic boundary conditions 

Under AP boundary conditions the energy of a charge configuration is again given by (2.38b). In this case, H. c ({Qr}) 
is as given by (5.6), but since now ir y = and ir x = ir, the energy of the charge configuration {Qj>} has an extra 
contribution n 2 JNM^ 1 from the first two terms in (2.38b), which is exactly the spin wave energy encountered in the 
discussion of the ferromagnet in Section 3. One cannot, however, without further inspection, identify this contribution 
as the ground state energy difference AE AP , since {Q%} is not necessarily the ground state any more. The remainder 
of the discussion depends on the value of £ x /£ y . We discuss first the case n£ x /£ y > 1. Curiously, in this case the 
mechanism to construct a lower lying state is the same as it was for R boundary conditions, viz. columnwise charge 
reversal. 

Reversing a single column leaves M~ 1 P X — but changes N~ 1 P y from to |, which annihilates the effect of n x = ir 
and therefore cancels the spin wave contribution w 2 JNM~ 1 . The same cancellation can be obtained by reversing the 
charges in an odd number of otherwise arbitrary columns. The energy cost is the excitation energy of the columns. It 
can be minimized by choosing an odd number of consecutive columns a, a + 1, ... a', such that the distances x a _i — x a \ 
and large as possible. Typically, they will again be of order £ x \og(M/£ x ). Hence we arrive for 

AE AP at the estimate 

Cy t X 

identical to AE R . This energy difference is less than the spin wave energy n 2 JNM^ 1 , and therefore the ground state 
will adjust to AP boundary conditions by a chiral excitation, and not by a spin wave excitation. In view of Eq.(1.3) 
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and our result for v c found above we now conclude that v s = v c = 1/(1 — tt£ x /£ v ). Hence at low temperature the spin 
and chiral correlation lengths, £ s and £ c , behave as 

6 ~ £c ~ r-V(i-^./<») (tt4/^ > 1) (5.11) 

In case < tt£ x /£ v < 1 (large spatial anisotropy, and, as seen above, a low temperature phase with long range chiral 
order), the situation is different. The minimum excitation energy of the columns is still (N/£ v )(M/£ x )~ 7 ' ex / e y , but 
this is higher than the spin wave contribution which it cancels. Therefore in this case the ground state {Qr} remains 
unchanged (and in the spin representation acquires only an additional global spin wave), and 

AE AP = n 2 JNM~ 1 (0 < ir£ x /£ v < 1) (5.12) 

In the thermodynamic limit M ~ N — > oo it follows (see Eqs. (1.2) and (1.3)) that v s = oo, which is most naturally 
(although not strictly necessarily) associated with a low temperature phase with power law decay of the spin-spin 
correlation. In this case the energy AE AP is not the result of a domain wall, but smeared out across the system. 

All of the above discussion is for an even number M' of charge carrying columns. We now comment on the differences 
that intervene when M' is odd. In that case the value of N~ 1 P y in (2.36) is ±i and is changed to (modulo 1) 
under replacement of P by AP boundary conditions. Hence the value of Qm/n does not change and (since M~ 1 P X = 
both before and after) we have AE AP = 0. The interpretation is that a global spin wave is "caught" in the system 
and cannot release its energy. When we replace P by R boundary conditions, the global spin wave terms disappear 
from the Hamiltonian and release their energy so that AE R is negative. The phenomenon of spin wave energy release 
was first observed by Kawamura and Tanemura [3] in their Monte Carlo simulations of an XY model with general 
disorder. It appears here analytically, but seemingly as a consequence of the distinction between even and odd M' 
and hence closely tied to the present example. But we shall see in the next section that it also appears in the general 
case. 



VI. GENERAL CASE 

On the basis of the experience gained we now attempt a theory, admittedly heuristic, for the general case. Let, as 
before, H c denote the standard lattice Coulomb Hamiltonian (2.30), and Ti R its counterpart (4.8) with R boundary 
conditions. Let in this section TL P denote the Hamiltinian H of Eq. (2.38b), and TL AP its counterpart with AP 
boundary conditions. We shall take the unknown ground state Q° of TL C as our reference state and denote its energy 
by Eq. Our strategy will be to try to determine how the ground states of H R and ii P ' AP differ from Q°. 

Since the remaining discussion concerns the low-temperature regime, it is convenient to work with charges ±| on 
the frustrated plaquettes and charges zero on all other plaquettes. This simplification, that was also made by Villain 
[12], leads to the problem of an Ising model on a random lattice (namely, the one composed of centers of the frustrated 
plaquettes), with logarithmically increasing antiferromagnetic interactions and zero magnetization. The logarithmic 
interaction will be screened in a way that we cannot precisely describe [it has to be if the energy per charge is to 
be finite] and it is reasonable to imagine that there exist some effective interaction decaying as a power law between 
spatially separated neutral sets of charges. 

Assuming that we know the ground state Q°, the next relevant question is what the lowest lying excitations 
are. Clearly all excitations can be described in terms of charge reversals with respect to the ground state Q , or 
alternatively, in terms of contours on the dual lattice (given any reasonable planar graph representing the nearest 
neighbor relations on the lattice of charges). In the spirit of Fisher and Huse [26] we shall consider the lowest lying 
excitations involving the reversal of order ~ L 2 charges and localized in an area of linear size ~ 2L around a preassigned 
point in space. Such excitations will be called droplets. Let their excitation energy scale as L~ x l Vc , where, since we 
are below d(, the exponent v c is positive. Due to the long-range forces, the droplets so defined need not constitute 
single domains. [For short-range forces they should and reduce to the droplets of Fisher and Huse [26] theory] A 
typical scale-!/ droplet may have to be represented by a set of disconnected contours. We shall assume, nevertheless, 
that it is compact enough so that there is always, typically, one main contour enclosing ~ L 2 charges. This amounts 
to assuming that a domain wall can be defined, whatever its width, around the reversed area. Excitations that would 
not fall in this category are, for example, those that consist of many small-size domains dispersed in the volume L 2 , 
or fractal excitations extending throughout the volume L 2 but having fractal dimension less than 2. Our hypothesis 
is that such reversals have excitation energies at least as high as those of the scale- i droplets. 

The ground states Q°' R and Q° of H R and 7i c , respectively, must differ at least by a contour going around the 
torus in the y direction, that is, by a scale- N droplet. For the energy difference between the two ground states we 
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therefore have 



E^-Eq-cN- 1 ^" (6.1) 

where c is a (positive or negative) random constant, and v c is as before the exponent of the chirality-chirality correlation 
function, see Eq. (1.2). This gives a way to determine this exponent. 

We now turn to the boundary conditions P and AP. Let Q a,p (Q°' AP ) be the ground state of H p (H AP ) and let 
E p (E AP ) be its energy. Since H P ' AP differ from H c by the addition of two quadratic terms, we have necessarily 
Eq' AP — Eq > 0. In the ground state Q° these additional terms take a random positive value of order 1 ( at least, if 
we assume that M~ 1 P X and N~ 1 P y do so in the ground state). This value can be reduced by reversing domains of 
charges in Q°. The ground state Q°> p (Q°' AP ) is now determined by a compromise between the minimization of H c 
and the minimization of the additional terms, which drive the total electric dipole moment to a specific value. Since 
this is a global biasing force, and since the excitation energies go down with increasing length scale, we expect that 
the compromise leads to a Q°' P (Q°' AP ) that differs from Q° by an excitation on the scale TV of the system. This 
excitation might not be of the same type as the droplet excitation discussed above, and therefore we put 

E P ' AP - E ~ c p ' AP N- 1 /"' (6.2) 

where c p and c AP are positive random constants and the exponent v s is also positive, it is the spin correlation length 
exponent of Eq. (1.2). But since, by hypothesis, the droplet excitations are the lowest lying ones that exist at each 
given scale, the excitation that leads to (6.2) is, at best, also a scale- N droplet excitation in which case v s = u c , or is 
possibly a combination of droplet excitations on smaller scales in which case one might have vj 1 < v~ x . Hence 

v- 1 < v-\ (6.3) 

Upon combining these considerations we arrive at the conclusion 

AE AP ~ (c AP - c p ) N- 1 /"* 

AE R ~ -c p N- 1/u ° - cN- 1/v " (6.4) 

in which c p,AP > and c may be of either sign with c = 0. It follows that AE R ~ —c p N~ 1 / Us , which is negative. 
This is the phenomenon of the release of spin wave energy as one passes from P to R boundary conditions. One sees 
furthermore that 

jAE^Pf 1 ' 2 ~ {AE R y 1/2 ~ N- 1 '^ (6.5) 

i.e. both energy differences scale with the same power of N . 

We now discuss the two possibilities implied by Eq. (6.3). If it holds with the equality sign, we are in the case of 
a single correlation length for the spin and the chiral variables, as found in the example studied. This case needs no 
further comment. If Eq. (6.3) were to hold as a strict inequality, then there are two different exponents v s and v c . 
This would imply a longer correlation length for the spin variables than for the chiral variables. Moreover, AE R is 
then not the appropriate quantity to determine v c . In fact, in this case v c cannot be found from a comparison of P,AP 
and R boundary conditions, but should result from a comparison of TL C and H R according to Eq. (6.1). Kawamura 
and Tanemura [3] extracted the exponent v s from a numerical determination of AE AP and v c from a postulated 
expression involving both AE AP and AE R . In this work there is a strict inequality that goes in the sens opposite to 
(6.3). 

The conclusion is that we find no evidence for chiral order extending on a longer length scale than spin order. 

VII. FURTHER COMMENTS AND CONCLUSIONS 

We have considered the two-dimensional XY spin glass with ±J interactions. In sections 2 through 4 general 
formulas are presented that exhibit the interplay between chiral and spin variables in determining domain wall 
energies on finite M x N lattices with various types of boundary conditions. In Section 5 we have considered, as an 
example, a specific type of disorder with infinite ranged correlations in the y direction. The spatial distribution of 
the frustrated plaquettes depends on two parameters £ x and £ y , whose ratio is an anisotropy parameter. The scaling 
properties of the domain wall energies with system size lead us to the following conclusions. 
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(i) XY spin glasses with critical temperature T = 

For £ x , £ y — > oo and for Tr£ x /£ y > 1, the system of Sec. 5 has only a T = critical point at which the chiral and 
spin correlation lengths diverge with the same exponent v s = v c = 1/(1 — tt£ x /£ v ). We add here without proof that 
the phenomenon of a T = transition with ^ s = v c holds in the entire region of the £ x £ y plane determined by 

44(4 -1) sin 2 ^- >1 (4>1 and 4 = 1,2,...) (7.1) 

This phenomenon is analogous to what we found on the one-dimensional ladder [6] and tube [25] lattices. It is 
distinct, however, from the scenario proposed by Kawamura and Tanemura [2, 3] and by Ray and Moore [5] for the 
two-dimensional uncorrelated random ± J XY spin glass, according to which one would have < v~ < v~ x . In order 
to rule out the possibility that our conclusions are restricted to a specific class of correlated disorder, we consider in 
Section 6 the general case of random disorder. The heuristic theory presented yields a different inequality, namely 
< v~ x < v~ x . We therefore find no evidence for chiral order extending on a longer scale than spin order in two 
dimensions. 

One may now go one step further and speculate that since we had v s = v c for d = 1, the mechanism uncovered 
above, which enables the ground state to accommodate to a change from P to AP boundary conditions by a low 
energy chiral excitation, is general for uncorrelated random ±J XY model, and that v s = v c for all d < de- 

(ii) random XY models with nonzero critical temperature 

This class encompasses the model of Section 5 when Eq. (7.1) is not satisfied. In this case there is a low temperature 
phase with long range chiral order and, very probably, power law decay of the spin-spin correlations. This is the same 
scenario that is believed to hold for a fully frustrated two-dimensional XY model [20-22], which is, as a matter of 
fact, recovered in the special limit £ x = £ y = 1, and which is characterized by v~ x < = v~ x . It cannot, therefore, 
be identified with the scenario v~ x < < v~ x that Kawamura and Tanemura [3] propose for the t/iree-dimensional 
random ±J XY model. 

Finally we comment on the transition line 44(4 — 1) sin 2 (7r/24) = 1 between the regions (i) and (ii). Clearly the 
distinction between these two types of ground state behavior is not due to the density of frustrated plaquettes, which 
is 1/44' but t° the nature of correlations in their spatial arrangement. The different behaviors can be distinguished 
by the ground state response to a change between P and AP boundary conditions. In region (ii), the ground state is 
clastic (it deforms continuously) , whereas in region (ii) it is "brittle" (it responds by the formation of a chiral domain 
wall). This transition appears in this work somewhat as a byproduct of the analysis, but it is interesting in itself and 
merits a separate study. 
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